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Abstract 



We compute the dispersion curves for neutrinos propagating in an extremely 
dense electroweak plasma, in the presence of very strong magnetic fields of 
order B < MlL/e. The neutrino self-energy is calculated in the one-loop 
approximation. We consider only contributions of the first Landau level to 
the propagator of the PL-bosons, and distinguish between motion parallel 
or perpendicular to the external magnetic field. We find that the neutrino 
dispersion curve for parallel propagation to the field suggests a superfluid 
behavior. An interesting analogy with fractional QHE is pointed out. We 
obtain a neutrino effective mass which increases with the magnetic field. 



I. INTRODUCTION 

The present paper is addressed to investigate the dispersion equation of massless neutral 
fermions, interacting with charged fermions and massive vector bosons, propagating in a 
medium at finite temperature and density, and in presence of an extremely strong magnetic 
field. 

The propagation of neutrinos in an electroweak plasma has been studied and the dis- 
persion equation for the quasiparticles was obtained ^HU- The spectrum found exhibits, in 
some extreme conditions, a superfluid behavior. 

In the present paper we consider the role of extremely strong magnetic fields as a possible 
mechanism for generating an effective neutrino mass in a very dense medium. We find again a 
superfluid behavior for the neutrinos moving parallel to the external magnetic field, provided 
it is strong enough. 

The propagation of neutrinos in magnetized media, assuming no dependence of the W- 
propagator on the magnetic field, has been done by computing two types of diagrams (bubble 
and tadpole) . We work with the effective action, the generating functional of irreducible 



one-particle inverse Green's functions. Thus, in the evaluation of the inverse Green's func- 
tion, whose zeros give the dispersion equation, the tadpole diagrams do not appear. The 
tadpole diagram is reducible and does not contribute to it. 
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For strong enough magnetic fields, a gas of charged bosons undergoes Bose-Einstein 
condensation ||. This suggests that it is the ground state of the bosons which play the main 
role. In our calculation we will take into account only the part of the W propagator which 
contains the ground state energy 0]. When the momentum is small and the magnetic field 



is high enough eB < M^r, so that the term — eB dominates, the main contribution 

to the propagator comes from the low momentum gauge bosons (IV-condensate) 0,0]. 

The expressions we derive for the neutrino dispersion equation are valid for strong fields, 
close to the limiting value M^/e. Such fields are at present conjectured to exist in the cores 
of neutron stars , and they may have also existed in the early universe, in which case the 
observed galactic and intergalactic magnetic fields are viewed as relics of huge primordial 



fields [13]- |T7]. However, on a more general basis, our results apply to any massless fermions 
interacting with vector bosons of mass M in strong magnetic fields eB < M 2 . 

The first part of this paper (section 2) presents some general expressions of Green func- 
tions in a magnetized medium, which will be used subsequently. In section 3, we find the 
mass operator in the limit of high magnetic field, after performing the sum over p 4 . These 
results are very general and apply equally well to the limiting cases of high temperature 
and high density. In section 4, we analyze the limit T — > 0, corresponding to a degenerate 
fermion gas (large fermion density). We also consider the role of VF-boson condensation. In 
section 5, the dispersion equations are discussed. We analyze two cases of neutrino prop- 
agation, parallel and perpendicular to the magnetic field. For parallel motion, an effective 
mass results only in the sense of lim^^o^ 7^ 0- When the motion is perpendicular to the 
magnetic field, a minimum appears in the dispersion equation at non-zero momentum and it 
is possible to define an effective mass in the more traditional sense of being a local minimum 
of the dispersion equation. 

Our results depend linearly on the magnetic field and, in the degenerate case, also linearly 
on the chemical potential fi. Other studies of the mass operator in QED in presence of a 
magnetic field || analyze different limits: high temperature, high density, and high magnetic 
field. The first two show similar dispersion equations || if one substitutes M = e 2 T 2 /8 with 
M = eV/(8vr 2 ). 



II. NEUTRINO SELF-ENERGY: GENERAL EXPRESSIONS 

The neutrino two-point inverse Green function in presence of a magnetic field reads as 

s-\k*) = p R (-i^k; + ^ w {k*))p L . (i) 

where k* x = k\ — ifi u 5^\, and Pr = (1 + 7s)/2 and Pl = (1 — 7s)/2 are the left- and 
right-handed projection operators. We recall QXj| that the following relation holds among 
the chemical potentials \i v = \i e — Hw~- Note that, in equation (P, the self-energy of the 
Z-boson is not included. It can be neglected since it does not interact with the magnetic 
field (at the one-loop level); it is of order g /2 /Mf , whereas the W term is of order g 2 eB jM^. 
The expression for the self-energy Tj W is the following, in configuration space: 

Z w (x,x') = -i-f^^ G e (x,x')D™(x - x') 7w , (2) 
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It represents the self-energy due to electron-W polarization. 

In Euclidean space and in the gauge = (0, Bx, 0, 0), the propagator of the electron is 



{{{ip4 - fj)nf4 + ipzl-z ~ m e )(a + ij n tp n + 0-_^„-iV>n-i) 



+ !/2V2aBn[7 + '0nV'n-i - T-V'n-i^]} (3) 



■exp[ipl(x 4 - x' A ) + ip 3 (x 3 - x' 3 ) + ip 2 {x 2 - x 2 )], 



where £ = VeB(xi + x Q ), £' = VeB(x 1 + x ), x Q = |§, a ± = l/2[l±a z ), j± = l/2[j 1 ±ij 2 ), 
a 3 = i/2[ji, 72], and p\=p\- i/iw^x- 

The W-propagator in a magnetic field has the form 



•'0n(OV'n(O ex P[*P4( a; 4 ~ Z 4 ) + W^2, ~ X3) + iPlfa ~ 4)]- 

where R± = \pf + E%, ± 2eB}~\ R° = [pf + E^,)" 1 , with E 2 v = M 2 v +p 2 3 + 2eB{n + 1/2); 

= 4^ = ^-srGS and 4^ = (G 02 is the fidd tensor of the SU (2)xU(l) 

electromagnetic external field). Concerning the gauge fixing term, we are taking D^ v in a 
transverse gauge which is expected to guarantee the gauge independence of the neutrino 
spectrum. 

The poles of D^ v are located at 



which is the ground state energy, and at 



= sjpl + w? w +2eB{n+ 1 -), (6) 

where n = 0, 1, 2, . . . with degeneracy (3 n = 3 — 5 0n . The ground state energy (H) is unstable 
for pg < eB — M^. The analog of the Euler-Heisenberg vacuum energy due to vector boson 
polarization is || 

1 r oo a+ ^7 
U w = -— — / -^e-^ieB csch(e£t)(l + 2cosh2e£t) - - - ~e 2 B% (7) 

Convergence of this expression is only possible for eB < Myy, i.e. the vacuum becomes 
unstable for eB > M^. This problem has been the subject of investigation mainly by 
Nielsen, Olesen and Ambjorn ||10|| , |fLl|| . In the last reference, a static magnetic solution of 
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classical electroweak equations, corresponding to a vacuum condensate of W and Z bosons, 
is found. It is valid above the critical value B c = / e. The vacuum bears the properties 
of a ferromagnet or an ant i- screening superconductor. 

We are interested in the Fourier transform of (||]). It requires rather long calculations 
involving the Fourier transform for two Hermite functions, which lead to functions of gen- 
eralized Laguerre polynomials. Eventually we find 

^ W ( k ) = |^(E / T^^G* + *s,rf + kl,n')J^)P L , (8) 
pi \ i 

where 

G e {P3 + h,Pi + h,ri) = ((p 3 + hf + (p 4 + k 4 ) 2 + m 2 e + 2eBn'y l 
and T, a /3 is a 4 x 4 matrix whose elements are the following (£12 = £21 = £34 = £43 = 0): 
£11 = #i(z(p4 + kl))T*,_ hn T n >- 1:n + Bx(-i(pl + kt))T*,_ hn T n ^ 1}7l 

+2A 1 (i(pl + kl))7Z, tn T n ,, n , 

£13 = -2Ai(p 3 + k 3 )T*, n T n/>n , 

£14 = -2iB 1 V 2eBnT*,_ ln T n ^ n - 2iCV 2eBnT*, n T n ,_^ n) 
£22 = B 1 (i(pl + kl))T*, iV T*, )n + Bti-iipt + fcl))T n *, in T n> 
+2A 1 (i(pt + kl))T:,_ 1>n T n ,_ hn , 

£ 23 = -2iB^2eBnT* n , J n ^ n - 2iB x v '2eBnT* n , _ x>n T n ., n - 2iC^2eBnT*,_^ n T n , >n , 

£ 24 = 2A 1 (p 3 + fc 3 )T n *,_ ln T n /_i >n , 
£31 = 2A 1 (p 3 + /c 3 )T^ n T n> , 

£ 32 = 2iB^2eBnT* n ,^ n T n ^ n + 2iC^2eBnT*,_ ln T n , , n , 

£33 = + /c*))T n *,_ ln T n /_i jn + £i(-i(>4 + A;4))T n *,„ lin T n /_ lin 

2A 1 (-z( P : + A;:))T* jn r ri , )n , 

£ 42 = 2Ai(p 3 + fc 3 )T n *,_ ln T„/_ lin , 

£ 44 = 5i(z(p: + fcI))T*, jn T n> + S^-z^ + kl))T:^ n T n ,, n 

-2A 1 (-z( P : + A;:))T n *,_ lin T n ,_ lin , 
where = + i?~)/2 + 2i2°, Si = C = -i/2(R + - iT), 

and £™ _n are the Laguerre polynomials. 
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III. NEUTRINO SELF-ENERGY IN THE HIGH MAGNETIC FIELD LIMIT 



Let us now consider the limit of an extremely strong magnetic field (eB close to, but 
smaller than, Mjy). The W propagator is dominated by the Landau ground state term 
(n = 0), which means we keep in @ only terms proportional to R~ . 

Furthermore, the condition eB > \j? e — m 2 implies that the only electron state which 
contributes to the mass operator is n' = 0. The sum J2^=o can be approximated to X^ =0 , 
where n M is the integer part of (p, 2 , — m 2 )/2eB, which is zero whenever p? e < eB, i.e., for 
most cases of interest. Hence, we concentrate on the case when both electron and W-boson 
are in the Landau ground state, and their quantum numbers are n = n f = 0. 

Keeping in mind the above approximations, equation (|8]) becomes 

^ W ( k ) = / dp 3 G e (p 3 + h,p 4 + h)^ aP )P L , (9) 



(2*) S „, 
J R-(p 3 ,P4)e- fe i/ eB S^, and R; 



E' 



a3 



(10) 



with E a/3 = R-(p 3 ,p A )e- k y e »Z' a p, and R~ = [M& + pf + p\ - 2eB)'\ G e (p 3 
\{p* A + kl) 2 + (p 3 + k 3 ) 2 + m 2 } 1 . The matrix E^o simplifies and takes the form 

((i(p% + k%) (p 3 + k 3 ) 0\ 



-(p3 + h) -i{p$ + kl) 

\0 0/ 

After performing in (|9|) the sum over p\ and taking the analytic continuation [k\ — > ik Q ) 
we get a function of the new variable k Q — p v . The singularities in this variable lead to the 
gauge-invariant, physically relevant spectrum. We have 

-z(&4 - iii v ) = (k Q — \i v ) = uj + iT 

where u is the energy and T the inverse lifetime of the neutrino quasiparticles. 
We obtain the following expression for E w : 

g 2 eB 



—[ u I 2 - h + I& 

^33 — ^11 



Mr 



Ml' 



g 2 eB 



[-k 3 I 2 + h- I 



5J- 



13 31 

where the integrals Ii can be written as (neglecting the vacuum contributions), 



(11) 



h = l% ((J °° ~ 2P3k ^ Ue ~ ^ ~ 2ujE ^e + Tip)) e- k±2 ^ B 



dp 3 



2EwQ 
dp 
2Q 



{J' OQ - 2p 3 k 3 )(n w - + n w +) + 2uE w (n 



w 



n 



w. 



h = I 2%> ( (J °° + 2p ^ nw ~ ~ nw+ ^ + 2K 



-k± 2 /2eB 



-k± 2 /2eB 



h = 

P3dp 3 



Pzdp 3 



{{J oo - 2p 3 k 3 )(n e + n f 



2EyyQ 



2E e Q 

7 ((J'oo + 2p 3 k 3 )(n w - + n w +) + 2E w u(n^ 



wu(n w + r% ^ 
2ujE e (n e -n p ))e- k±2/2eB 

-k± 2 /2eB 



n 



w. 



(12) 
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with 



Q = {J 00 - 2p 3 k 3 ) 2 - Au 2 E 2 e 
Q' = {J' 00 + 2p 3 k 3 ) 2 -Auj 2 E 2 w 



and 



Joo = z 1 - eB - m 2 + Mw, 
J' oa = Zl + eB + m 2 e - Myy- 



In the above formulae E e = J ((p 3 + k 3 ) 2 + m 2 ) and Ew = y(p 2 + M^ — eB) are the 
Landau ground state energies of the electron and the I^-boson, respectively, whereas 



are respectively the distribution functions of the electrons, positrons, W~ and W + in our 
plasma. 

The mass operator given by the expressions ® is general and holds also in the high 
temperature and high density limits. The branch points in the denominators Q and Q' can 
be identified as thresholds for neutrino absorption in the plasma. We postpone a careful 
study of the analytic properties of the mass operator and their implications for the dispersion 
equation to future work. 



In this section we consider the case of degenerate electrons (formally equivalent to the 
limit T — > 0). Our results are of interest in the theory of neutron stars as well as in the 
early universe, when there might have been magnetic fields eB ~ (10 22 gauss). In the 
degenerate case, the distribution of the electrons is just a step function, n e = 9(fi e — E e ), 
and there are no positrons left, n p = 0. Charge neutrality is ensured thanks to some 
W + background. From the behavior of the distributions of W^ 1 0|| it has been argued 
that iy-condensation in the presence of a magnetic field may indeed take place. Thus, 
the distribution of W + can be approximated by (2Tr 2 )8(k 3 )^- (Nw is the total density of 
VT-particles in the medium), and for the excited states nw+ = 0. 

In this limit, equations become 



n 




-l 



n w - tW + = [ e ^=FM^)/3 _ i] 



l 



IV. DEGENERATE CASE 





2(2vr) 2 
< 2 N w e- hA? f' 



(2E W ) k 3 + 2p 3 k 3 + m 2 e -uo 2 - 2uE w - d? ' 




2(2tt) 2 
2 eBe~ k±2/2eB 




(13) 



where d 
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V. DISPERSION EQUATION 



Before solving the dispersion equation, note that we work far from the thresholds for 
neutrino absorption. In order to get the dispersion equation we must solve 



det(-z 7 A + S vv ) =0 



which yields 



(k 2 3 - u 2 ) 



E 33 )o; 



2h{h - £13) + kl - 2uo 



^11^33 + (&3 
UjCEii — E33) 



s 13 ) 2 
= 0, 



(14) 

(15) 
(16) 



Let us remark that the limit eB — > does not mean that Tj W = 0. The dispersion 
equation (|16D, when k 3 — > k± — > 0, leads to a value for u different from zero; this value 
corresponds to a sort of "effective mass" |]J proportional to eBfi e /d 2 . This means that, for 
fixed electron density, if the magnetic field grows up to near M^/e the "effective mass" grows 
too. As pointed out before, the vacuum is unstable for eB > Myy We shall show below 
that when the motion is perpendicular to the magnetic field, the "effective mass" becomes 
a mass in strict sense, since it is also a minimum of the dispersion curve: in formulas, 
8lu / dk±\ k±=0 = and 8 2 lu /dkj_\ k±=0 > 0. 

Since we are considering the degenerate case, the contribution of terms containing the 
W-bosoia distribution function can be safely neglected: for huge magnetic fields, Nw ~ 
C « eB, (C is the W-condensate), whence only Ii and I4 contribute in (|16|). 

In order to solve numerically equation (H)), we distinguish two cases: motion parallel to 
the field (k± — > 0) and motion perpendicular to it (k 3 — > 0). This equation involves only Ii 
and 1 4: 



(k 2 3 -LU 2 + k 2 ± ) 2 - (k 2 3 -LU 2 + ki)(2k 3 h) + (kl - LU 2 )lf = 0, 



(17) 



A. Motion parallel to magnetic field 

Equation ([IT]) for neutrino propagation along the magnetic field becomes 

(kl - u 2 ) 2 + (k\ - u 2 ) [21* - 2k 3 h] = 0. (18) 

Figure (1) shows the neutrino dispersion curves in this case, having fixed \x = 100 m e and 
eB = 0.9M^. It has two branches. One of them corresponds to the light cone. The second 
one arises due to the magnetic field and the finite density. The non-zero intercept at k 3 = 
is a sort of "effective mass" (for our typical values it is approximately 0.92 m e ). Besides, 
the curve has a gap (0.68 m e ). 

Interestingly, the curve shows a close analogy to the collective excitations arising in the 
fractional quantized Hall effect. In paper a theory of the excitation spectrum in the 



fractional Hall effect analogous to Feynman's theory for the excitation spectrum of superfluid 
helium was proposed. A magneto-roton minimum for the collective excitation spectrum was 
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found, which has a remarkable analogy with the minimum obtained in the present case, 
when the neutrino propagates parallel to the magnetic field. 

It is possible to interpret the gap of the quasiparticle spectrum as the symptom of a 
superfluid behavior. A similar interpretation has been done in the case of the dispersion 
of neutrinos in a hot medium without magnetic field. But here we must observe that 
the neutrinos interacting with the W-s and electrons must align their spins also along the 
magnetic field, leading to weak coupling in pairs, and to condense. 



B. Motion perpendicular to the magnetic field 



When the neutrino moves perpendicularly to the field, we get the following expression 
for the dispersion equation 

-uj 2 (-uj 2 + 2I 2 ) + k\{k\ -2c; 2 ) = 0. (19) 

The numerical solution also yields two branches (figure 2); one of them gives the same 
"effective mass" as in the parallel case. This is a mass in the proper sense since it is the 
minimum of the dispersion equation. 

In spite of the same value for the effective mass, the dispersion curves at zero momentum 
have different slopes for motion parallel and normal to the field. The behavior of both curves 
is quite different in both cases. This conclusion is to be expected since the magnetic field 
produces an anisotropy in the system and the motion in these two directions have different 
physical properties. However, the most notable result here is the behavior of the effective 
mass ra v e fj ~ A i e e-B/(M^ / — eB), which increases without bound as eB — ► M^. For fields 
eB > M$y, the neutrino magnetic mass problem, requires further research along with the 
Higgs mechanism in external fields, taking into account the results of refs. [pi], [11 . 
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VII. APPENDIX 



We present here the result of the calculation of Ii far from the thresholds. Taking into 
account the condition 

k\ — iJ 1 + d 2 > /J e (k 3 — to), 

we get 



S 



g 2 r 1 
L = 2{2n) 2 kj - uu 2 + d 2 ' 
9 2 N W 1 
2(2tt) 2 d{kl - (to + d) 2 ) ' 

g^V 1 

2(2tt) 2 fc§ - (w + d) 2 ' 
h = h 
h = 0. 
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Fig 1. Curve of propagation of neutrino paralell to the magnetic field 




Fig 2. Curve of propagation of neutrino normal to the magnetic field 
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